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Abstract: 

Using the method of point canonical transformation, we derive some exactly solvable 
rationally extended quantum Hamiltonians which are non-Hermitian in nature and whose 
CN ■ bound state wave functions are associated with Laguerre or Jacobi-type Xi exceptional 

orthogonal polynomials. These Hamiltonians are shown, with the help of imaginary shift of 
co-ordinate: e~ ap xe ap — x + ia, to be both quasi and pseudo-Hermitian. It turns out that 
the corresponding energy spectra is entirely real. 
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INTRODUCTION 



Since the discovery [l|, 0] of X\ exceptional orthogonal polynomials (here after EOPs) in math- 
ematical physics there has been renewed interest in the analysis of exactly solvable shape invariant 
quantum systems. Unlike the classical orthogonal polynomials, these new polynomials have the 
remarkable properties 0] that they still form complete sets with respect to some positive definite 
measure, although they start with degree n > 1 polynomials instead of a constant. Laguerre and 
I/-) | Jacobi-type X\ EOPs have made their appearance in the bound state wave functions of the quantum 

systems with both constant 0, 0] and position-dependent mass [f|. These quantum systems are 
shown [7j, with the help of reducible second order supersymmetric transformation, to be rationally 
extended version of conventional ones associated with the classical orthogonal polynomials. This 
supersymmetric transformation also explains the isospectrality of the conventional and rationally 
extended potentials. Subsequently, EOPs are generalized to higher co-dimension @-10] and to 



multi- indexed systems 11, [lj] and associated shape invariant Hamiltonians are reported 13- 17]. 
Some properties of these polynomials are studied in ref. [1J]. EOPs are also used in connection with 
discrete quantum mechanics 0, Dirac and Fokker-Planck equations 21], pre-potential approach 
[22I ]. information entropy 23]. quantum Hamilton Jacobi formalism (24)], dynamical breaking of higher 



order supersymmetry [25[ and quasi-exactly solvable problems [26f] . However, the application of these 
new polynomials to the non-Hermitian quantum systems is not reported so far. 

Non-Hermitian Hamiltonians are important due to the fact that, despite being non-Hermitian in 



nature, these operators may constitute unitary quantum mechanical systems [2714291] . Non-Hermitian 
parity-time {VT) symmetric Hamiltonians possess real discrete energy eigenvalues if the corresponding 
eigenfunctions are also VT symmetric, otherwise the eigenvalues occur in complex conjugate pairs 
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3C| • "PT-symmetric Hamiltonian having all eigenvalues real is connected to the existence of a positive 



definite inner product (t^,rj<p) which render the Hamiltonian H to be pseudo-Hermitian [311434 1 
H< = rjHr] -1 , where the Hermitian linear automorphism r\ : T~L — > T~L is bounded and positive definite. 
Another equivalent condition for the reality of the energy spectrum H is the quasi-Hermiticity 
28|, 35-37], i.e. the existence of a invertible operator p such that h = pHp~ l is Hermitian with respect 
to usual inner product {ip\4>)- Quasi- Hermitian Hamiltonian shares the same energy spectrum of the 
equivalent Hermitian Hamiltonian h and the wave functions are obtained by operating p -1 on those 
of h. Most of the analytically solvable non-Hermitian Hamiltonians are constructed by making the 
coupling constant of the known exactly solvable potentials imaginary 3^-41]. In some other cases 



[42j-|44j], the coordinate is shifted with an imaginary constant. Several of these classes of Hamiltonians 
are argued to be pseudo-Hermitian under rj = e~ ap 45, Hf]]. For a real a and p = —i-^, the operator 
rj shifts the coordinate x to x + ia. 

The goal of this letter is to generate some rationally extended Hamiltonians which are non- 
Hermitian in nature and whose bound state solutions are associated with Laguerre or Jacobi-type 
X\ exceptional orthogonal polynomials. By 'rationally extended Hamiltonians' we mean those which 
are the extensions of the well known Hamiltonians by addition of some rational functions. The 



method of point canonical transformation (PCT) [47J, |48(], which consists of transformation of the 
initial Schrodinger equation to a differential equation of some special function, has been used here 
to achieve our goal. The non-Hermiticity enters, in a natural way, into the potentials through 
the purely imaginary constant of integration appears in PCT. We also show, with the help of a 
similarity transformation, that the new non-Hermitian Hamiltonians obtained here are quasi as well 
as pseudo-Hermitian. In particular, it has been identified that the positive definite operators p = e~? p 
and rj = p 2 play the roles of a quasi and pseudo-Hermitian operators respectively. 



II. QUASI-HERMITIAN HAMILTONIANS ASSOCIATED WITH LAGUERRE OR JACOBI 

TYPE Xt EOPS 

Here we use the method of PCT to derive some exactly solvable non-Hermitian Hamiltonians whose 
bound state wave functions are associated with Laguerre or Jacobi type X\ exceptional orthogonal 
polynomials. For this we first briefly recall the method of point canonical transformation. 



In PCT approach [471. l48j . the general solution of the Schrodinger equation (with h = 2m = 1) 



HiP(x) = -^rP + V(x)iJj{x) = Ei){x) (1) 

can be assumed as 

V(z) ~ f{x)F(g{x)) (2) 
where F(g) satisfies the second order linear differential equation of a special function 

d 2 F dF 

^ Y + Q{g) ir +R{g)F(g)=Q. (3) 

Substituting the assumed solution ip(x) in equation ([I]) and comparing the resulting equation with the 
equation (J3j) one obtains the following two equations for Q(g(x)) and R(g(x)) 

q" 2 f 

Q{9) = ^ + JJ (4a) 
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respectively. After some algebraic manipulations above two equations reduces to 

/(sWOO" 1 / 2 e i/0W *, (5a) 

9"' 3 fg"\ 2 , /2 / ldQ 1 2 



Now, we are in a position to choose the special function F(g) (consequently Q(g) and R(g)). The 
equation (|5bp becomes meaningful for a proper choice of <?(a;) ensuring the presence of a constant term 
in the right-hand side which connects the energy in the left-hand side. The remaining part of equation 
(|5bp gives the potential. Corresponding bound state wave functions involving the special function F(g) 
are obtained with the help of equations §Z§ and (|5a|) . as 

^x)~g'(xy^e hlQ{9)d9 F(g(x)). (6) 

Here, we choose the special function to be the exceptional X\ Laguerre polynomial viz, F(g) oc Ln\x). 
For real a > and n = 1,2, 3..., these polynomials Ln\x) satisfy the differential equation [l[ 

d 2 y (a-q)(a + a + l) dy ( x-a \ 

x ~r^ t + hn-1 h/ = 0. (7) 

qar x + a ax \ x + q J 

The polynomial L^(x) has one zero in (— oo, —a), remaining n — 1 zeros lie in (0, oo) . Moreover, these 
polynomials are orthonormal [2] with respect to the rational weight W = 



/„ (i + a)2 » * ' v ' (n-1)! 
The expressions for Q{g) and -R(g), corresponding to the choice F(g) = Ln (x), are given by 

ni \ (g-a)(g + Q + l) „, v g-a , n-1 

= 7 — ; — \ , = -7 — : — T + • ( 9 

Using them in equation ()5bp . we have the expression for — V(x) as 

F_W y 9"' 3/VY (2nq + a 2 - q + 2 )g' 2 g ' 2 (a 2 - l)g' 2 2g' 2 g« 

KW 2 5 ' 4 vS 'J 2a<? q( 5 + a) 4<? 2 ( 5 + q) 2 4 1 Uj 

At this point we choose g' 2 / g = constant = k 2 , k G E — {0}, which is satisfied by 

g( X ) = - A (kx + df, (11) 

where d is an arbitrary constant of integration. Here two cases may arise, namely, d = and d ^ 0. 
Without loss of generality we can choose, for the moment, d = 0. For this choice, substituting g(x) in 
equation (fTOjl and separating out the potential and the energy, we have 

V(r\ - ^<r 2 -I- (a 2 - l \ 1 4- „ 4fc2 ,32afc 2 

V {X) — lg X f ^q 4j? + P?+S (fc 2 :r 2 +4a) 2 ' 

(12) 



fc 2 (2n+a-l) 



■>n — 2 
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The potential V(x) is singularity free in the interval < x < oo. The same potential has earlier been 
reported in ref.[4j. It has been shown that the potential V(x) is the extension of the standard radial 
oscillator by addition of last two rational terms. Such terms do not change the behavior of the potential 
for large values of x, while small values of x produce some drastic effect on the minima of the potential. 
The normalized wave functions corresponding to the potential can be determined, in terms of Laguerre 
X\ EOPs, using equations ([6]), ([8]) and 



(n - 1)1 k 2a+2 



2 2a-3 ( a + n -)Y(a + n - 1) J k 2 x 2 + 4a 



as 



/ I —\ • n- .1.2.3... 



(13) 



It is worth mentioning here that the choice d = in equation (jlip always gives rise to Hermitian 
potential. Nonzero real values of d do not make any significant difference in the potential and its 
solutions. The non-Hermiticity can be invoked into the potential only if d is purely imaginary. We set 
d = ie, e G K — {0}, and g(x) = \{kx + ie) 2 for which the potential reduces to 



k 2 (kx + ie) 2 k 2 (a 2 - |) 

V{X) = — h — . . Tn + 



Ak 2 



32ak 2 



(14) 



16 ' (kx + ie) 2 (kx + ie) 2 + 4a \{kx + ie) 2 + 4a] 2 
The above non-Hermitian potential is free from singularity through out the whole real x axis. Since the 
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FIG. 1: (a) Plot of the real (solid line) and imaginary (dashed line) parts of the quasi-Hermitian Potential 
V(x) associated with Xi Laguerre EOPs. (b). Plot of the corresponding equivalent Hermitian potential 
V(x) (thick line) and square of the absolute value of its lowest two wave functions. Here we have considered 
e = 1.2,a= 2,k= 1.75. 



energy E n has no dependence on d, the non-Hermitian Potential V{x) also shares the same real energy 
spectrum of V{x). This requires further explanation. In the following, we show that the potential V(x) 
is actually quasi-Hermitian. For this we define the operator 



d 

p = -i— 

ax 



which has the following properties 



pxp 



ie 



ppp 



P, Pf(x)p 1 = f 



(15) 



(16) 



In other words, the operator p has an effect of shifting the coordinate x to x — For the proof of the 
results (|16|) . readers are advised to follow the reference 45|]. For this operator we have the following 
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similarity transformation 



P V(x)p- 1 =v(x-j)= V(x). (17) 



This ensures that the non-Hermitian Hamiltonian corresponding to the potential V(x) is quasi- 
Hermitian. The equivalent Hermitian potential V(x), which corresponds to d = 0, is given in equation 
(fT2|) . It is very easy to show that the positive definite operator i] = p 2 satisfies ^(x)?? -1 = V^(x) 
ensuring the potential to be pseudo-Hermitian. The potential V(x) also satisfies V*(— x) = V{x) 
and hence is "PT-symmetric. The wave functions of the potential V{x) can be determined by 
4> n {x) = p~ l ip n (x) = tf) n (x + f ). 

In figure 1(a), we have shown the real and imaginary parts of the potential V(x) given in (|14|) . while 
figure 1(b) shows its equivalent Hermitian analogue V{x) given in (|12p . Using first two members of 
exceptional X\ Laguerre polynomials L^\x) = — x — a — l,L^\x) = x 2 — a(a + 2), we have also plotted 
in figure 1(b) the absolute value of first two wave functions given in (| 1 3 j) . 

Next we choose F{g) to be Jacobi-type X\ EOP, Pi a ' b \ which is defined for real a,b> —1, a ^ b 
and n = 1, 2, 3... In this case the expression for Q(g) and R(g) are given by [2| 

. (a + b + 2)g + a- b 2(6 - a) , (b - a)g - (n + a + b)(n - 1) (a - b) 2 

Q{9) = ; — ? 71 ^ 1 ' R \9> = — 



g 2 (b — a)g — 6 — a' 1 — g 2 (b — a)g — b — a 

(18) 

Using these expressions in (|5b|) and choosing g' 2 /(l — g 2 ) = constant = k 2 (k ^ 0), we have 

g{x) = sin(/cx + d). (19) 

Like the exceptional Laguerre polynomials, the choice d = gives rise to the potential, energies and 
corresponding bound state wave functions, as 



v 1 4 2 a+b-(b-a) smkx [a+b-(b-a) smkxY ' 

E n = ^(2n + a + b-l) 2 , 



(20) 



and 



b , 1 



, . , (1 — sin kx) 2 + 4 (1 + sin kx) 2 + 4 

Vv. x ~ ^ 77 x ■ , 4 a ' 6) skifcx , n= 1,2,3... 21 

a + b — (b — a) sm kx 

respectively. The above periodic potential V(x), which is free from singularity in the interval — ^ < 
x < t^, can be interpreted [1] as the rational extension of the standard trigonometric scarf potential 
which is associated with classical Jacobi polynomials. The wave functions in equation (|2ip are regular 



41( iff a,b> -1/2. 



Here, the non-Hermitian potential corresponding to the choice d = ie is obtained as 

) 

(22) 



V(x) = k ^ 2a ~^ 2b — — sec 2 (kx + ie) — k ^ 2 a ^ sec(kx + ie) tan(A:x + ie) 



2k 2 {a+b) 2k 2 [(a-b) 2 -4ab] 

a+b-(b-a) sin(kx+ie) i" [a+f>-(b-a) sm(kx+ie)} 2 

This potential V(x), which is defined on whole real line, also shares the same real eigenvalues of 
the potential given in (|2Up . Like the rationally extended radial oscillator the above non-Hermitian 
potential is also quasi-Hermitian under the the operator p defined in (|15p . The corresponding 
equivalent analogue is the one given in equation (|2Up which corresponds to the choice d = 0. The 
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(a) 



(b) 



FIG. 2: Plot of the real (solid line) and imaginary (dashed line) parts of the quasi-Hermitian Potential V(x) 
associated with X\ Jacobi EOPs. (b). Plot of the corresponding equivalent Hermitian potential V(x) (thick 
line) and square of the absolute value of its lowest two wave functions. Here we have considered a = 1.75, b = 
3, k = 1.25,e= 1. 



potential V(x) also fulfills the requirement of PT-symmetry i.e. V*(— x) = V(x), only if a = b. 
However, if we consider the other solution g(x) = cos(kx + d) of </ 2 /(l — g 2 ) = k 2 , the corresponding 
potential V(x) becomes PT-symmetric for all real values of a and b. The wave functions of V{x) 
can be determined by operating p^ 1 on ip n given in (|2ip . Here, we have not considered the complex 
values of a, b because this will give rise to the exceptional Jacobi polynomials with complex indices 
and complex arguments. The orthogonality properties for such complex polynomials may depend on 
the interplay between integration contour and parameter values. 

In figure 2(a) we have plotted the real and imaginary parts of the potential V{x). The corresponding 
equivalent Hermitian analogue V(x) and square of its first two wave functions are plotted in figure 

(a,b) _ _x _ 2+a+b 



2(b) We have used the expression of first two members of Jacobi type X\ EOPs, P{ = — | — 2 (a-b) 
and P 2 (a,fe) = -^±f±^x 2 - a2+ 2 { t- h ) +b) % ~ to P lot the S( l uare of the wave functions. 



III. SUMMARY 

In summary, we have generated some exactly solvable non-Hermitian Hamiltonians whose bound 
state wave functions are associated with Laguerre and Jacobi-type X\ exceptional orthogonal polyno- 
mials. The Hamiltonians are shown, with the help of imaginary shift of coordinate, to be both quasi and 
pseudo-Hermitian. The imaginary shift of the coordinate enables us to make the potentials singularity 
free throughout the whole real axis. The obtained potentials enlarge the class of analytically solvable 
non-Hermitian potentials. In addition, the non-Hermitian rationally extended trigonometric scarf po- 



tential might has potential application in "PT-symmetric optical lattice 49, 50]. It is to be noted here 
that the other choices of g(x) in the expression E — V{x) associated with Laguerre and Jacobi EOPs 
give rise to the several other exactly solvable Hermitian as well as quasi-Hermitian extended potentials. 
But in all these cases we have to redefine the parameters carefully so that n dependent term appears 
only in the constant energy. 



We emphasize that analogous study 5l[ can be made to the case of solvable Hamiltonians associated 



with exceptional orthogonal polynomials of higher co-dimension and multi-indexed polynomials. 
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